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1. Introduction 

Solitons are particle-like waves that propagate in dispersive or absorptive media without 
changing their pulse shapes and can survive after collisions. Various types of optical 
soliton phenomena have been studied extensively in the area of nonlinear optical physics. 
These include the nonlinear Schrodinger solitons in dispersive optical fibers, spatial and 
vortex solitons in photorefractive materials/waveguides, and cavity solitons in resonators 
P. Solitons in optical fibers were first predicted by Hasegawa and Tappert |2J in 1973 
and were first observed experimentally by Mollenauer et al. |3] in 1980. Since then, the 
idea of using solitons for long-haul optical transmission has been an attractive research 
area with rapid progress 4.. In 1986, Gordon and Haus pointed out that the spontaneous 
emission from optical amplifiers in the optical link can cause the timing jitter variance 
of the solitons to be proportional to the cubic of the propagation distance This 
effect, known as the Gordon-Haus effect, will place an upper limit on the achievable bit 
rate in the long-haul soliton communication systems, even though it was found out later 
that one can somehow overcome this limitation by using guiding filters U\ ■ 

The spontaneous emission effects on soliton propagation already calls for a rigorous 
quantum theory for soliton propagation, even though it may still be possible to treat 
them semi-classically. In 1987, Carter and Drummond et al. |H1 E] used the positive- 
P representation approach to transform the quantum nonlinear Schrodinger equation 
into stochastic nonlinear equations with noise terms. They then solved the stochastic 
nonlinear equations numerically and showed that solitons are quantum mechanically 
squeezed during propagation. The squeezing ratio of the quantum solitons was later 
calculated with the inclusion of the homodyne detection scheme [TO]. Since then, the 
quantum theory of traveling-wave optical solitons has been intensively developed and 
several approaches have been successfully carried out. The studied soliton systems 
include the family of nonlinear Schrodinger solitons [HI ^JJ El El E] as well as the 
self- induced-transparency solitons The quantum theory of nonlinear Schrodinger 
solitons has also been directly applied to study the soliton communication problems 

PS EH EE]. 

The quantum theory of nonlinear Schrodinger solitons turns out to be an very 
interesting one. It was sometimes less known that the quantum nonlinear Schrodinger 
equation is the evolution equation of one-dimensional bosonic system with 5-function 
interaction under the second quantization framework By solving the problem in 
the Schrodinger picture using the Bethe's ansatz method, one can construct bound-state 
eigensolutions which are closely related to the soliton phenomenon. In 1989, Lai and 
Haus noticed this correspondence and constructed the exact soliton states by the Bethe's 
ansatz method [T2]. At the same time they also developed an approximate nonlinear 
analysis based on the time-dependent Hartree approximation [TT| . In the subsequent 
year, a quantum soliton perturbation theory based on the linearization approach was 
successfully carried out by Haus and Lai ^2] ■ The squeezing ratio of the soliton parts can 
be predicted by a simple analytic formulation and the optimal local oscillator pulse shape 



Quantum Theory of Fiber Bragg Grating Solitons 



3 



can be determined. This quantum soliton perturbation theory was later re-formulated 
to make it more clear ^7] and eventually led to the development of the back-propagation 
method [2 3 , which is a general numerical method for calculating the quantum noises of 
nonlinear pulse propagation problems based on the linearization approximation. The 
aim of the present paper is to further pursue the theoretical development along this 
line by generalizing the theory of quantum solitons to bi-directional nonlinear pulse 
propagation problems and to study the quantum properties of Bragg grating solitons in 
particular. 

On the side of experimental progresses, many experiments for actually measuring 
the soliton squeezing have been carried out since 1990, with the advance of stable pulse 
lasers and high quantum efficiency detectors at the optical communication wavelengths. 
To provide a stable experimental setup against the environmental fluctuations, Shirasaki 
and Haus proposed to use a Sagnac loop in soliton squeezing experiments for measuring 
the quadrature squeezing [2T]. The scheme is actually a squeezed vacuum state 
generator that can output squeezed vacuum states. Soliton squeezing from a Sagnac 
fiber interferometer was first observed by Rosenbluh and Shelby in 1991 (1.7 dB below 
the shot noise) [22J. Using the same setup but at the different wavelength, Bergman 
and Haus also succeeded to observe squeezing with non-soliton pulses (4.9 dB below 
the shot noise) [23J in the same year. Since then, larger quadrature squeezing from 
fibers has been obtained with a gigahertz Erbium-doped fiber lasers to suppress the 
guided acoustic-wave Brillouin scattering (GAWBS) and a 6.1 dB noise reduction below 
the shot noise has been reported [23|. One can also produce amplitude squeezing by 
using solitons with an external spectral filter or with an imbalance Sagnac loop 
interferometer (SHj. Squeezing generation with an imbalance Sagnac loop interferometer 
in the normal- dispersion region has also been demonstrated j^Z]. Recently, quantum 
correlation between different spectral components of the squeezed soliton has also been 
studied extensively in the literature [2H|- The above brief description summarizes the 
exciting development that has been carried out for squeezing generation using optical 
fibers. 

In recent years, with the advance of new fabrication technologies, it becomes more 
feasible to actually utilize one- or higher dimensional periodic dielectric structures 
(or especially the photonic bandgap crystals) j2H] for modifying the properties of the 
spontaneous emission as well as the propagation of waves. The fiber Bragg gratings 
(FBGs) are one-dimensional periodic structures with weak index modulation. The 
simplest way to induce 1-D Bragg gratings inside optical fibers is with the side- 
illumination of the UV interference lights. The FBG formed this way can be viewed 
as a 1-D photonic bandgap crystal for the guiding mode of the single-mode fiber [30J. 
Interesting soliton phenomena known as the fiber Bragg grating solitons [HU E2| can be 
found if the fiber Bragg grating has the third-order Kerr nonlinearity. The propagation 
of optical pulses inside the nonlinear FBG can be described by the nonlinear coupled 
mode equations (NCMEs). Intuitively the fiber Bragg grating soliton is formed when 
the input pulse has the suitable pulse-width and peak intensity such that the nonlinear 
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Kerr effect is large enough to compensate the high anomalous dispersion near one of the 
bandedges of the FBG. 

From the theoretical point of view, fiber Bragg grating solitons belong to the class 
of bi-directional pulse propagation problems, where the quantum theory is still lack of 
enough consideration. Most of the previous studies on fiber Bragg grating solitons have 
been on the classical effects and there is almost no reported result on their quantum 
properties. It is the aim of this study to bridge this gap by developing a general quantum 
theory for bi-directional pulse propagation problems and particularly to apply the theory 
for studying the case of fiber Bragg grating solitons. It will be shown that the fiber Bragg 
grating soliton pulses will quantum-mechanically get amplitude-squeezed after passing 
through the fiber grating and the squeezing ratio will be calculated theoretically. We use 
the linearization approach to study the quantum effects of FBG solitons by extending the 
back-propagation method previously developed [20] to the case of nonlinear bi-directional 
propagation problems. By following the same spirit of the back-propagation method, 
we will first derive a set of linear adjoint equations from the linearized nonlinear coupled 
mode equations in such a way that any inner product between the solutions of the two 
equation sets is conserved during the time evolution. In this way, the variance of the 
measured operator as well as its squeezing ratio can be calculated readily for a given 
measurement characteristic function. The squeezing ratio of FBG solitons will be found 
to exhibit interesting relations with the fiber grating length as well as with the intensity 
of the input pulse. 

In contrast to uniform fiber gratings, nonuniform fiber gratings with chirp and/or 
apodization have shown some potentials for pulse compression applications - The 
theory and experiment for Bragg grating solitons propagated in apodized FBGs have 
also been developed and demonstrated [23 ESI- Some new results on slowing down 
FBG solitons in apodized fiber gratings have also been carried out 36J. In the present 
paper the quantum fluctuations of FBG solitons in nonuniform fiber gratings will also 
be studied. We find that one can compress the soliton pulsewidth and enhance the 
squeezing simultaneously by using an apodized grating, as long as the solitons evolve 
adiabatically. 

The paper is organized as follows: in section [2] the model of the nonlinear coupled 
mode equations for bi-directional wave propagation in a uniform FBG is reviewed. In 
section El we use the back propagation method to calculate the quantum fluctuations of 
FBG solitons and the characteristics of the squeezing ratio are discussed. In section 
0] the quantum fluctuations of FBG solitons in nonuniform fiber Bragg gratings are 
studied. Finally, a brief conclusion is given in sectional 
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2. Fiber Bragg Grating Solitons 

We start the derivation from considering the linear wave propagation problem in a one 
dimensional periodic structure: 

d 2 E n\z) d 2 E _ 

'dz 2 ~ c 2 dt 2 ~ ( ' 

Here the dielectric constant n 2 (z) = n 2 + e(z) is a periodic function of the propagation 
distance, with the spatial average refractive index n and the period A. We are interested 
in the light field at the frequency near the Bragg condition ujq = k c/n, where c 
is the light velocity in free space, and k = n/A. For FBGs, we can expand the 
periodic index perturbation function e(z) by the Fourier series and only keep the 
phase-matching ±1 order terms, i.e. e(z) = 2eo cos(2k z). One then decomposes 
the light field into the forward (U a ) and backward (U b ) propagation pulses, E(z,t) = 
U a (z,t)e~ l( - u;t ~ koz ) -\-Ub(z,t)e~ l ( ult+koz * > + c.c., and obtain the following linear coupled mode 
equations: 

1 d d 

— ftU a (z, t) + — U a = i5U a + mU b (2) 
X d d 

-Qj_U b (z, t) - —U b = i5U b + mU a (3) 

where v g = n/c is the group velocity of the pulses, 5 = u — u is the wavelength detuning 
parameter, and k = u e/2nc is the coupling coefficient. In Fig. ^one can see that the 
dispersion relation of this set of coupled-mode equations has a bandgap at the frequency 
(wavelength) that satisfies the Bragg condition k = ko- 

If the third order nonlinearity of the optical fiber needs to be taken into account, 
one can model the problem by using the above coupled mode equations that describe 
the coupling between the forward and the backward propagating waves in a uniform 
FBG with the addition of the self-phase and cross-phase modulation effects. With the 
nonlinear terms, this set of NCMEs has analytical soliton solutions for the case of infinite 
grating length, as is shown by Aceves and Wabnitz with the introduction of the massive 
Thirring model [SI]. However, for nonlinear FBGs of finite length, no analytic solution 
can be found. So in our studies we directly use the finite difference numerical simulation 
method with the parameters based on the first experiment reported in the literature 
|32"] . We consider a 60 ps full width at half maximum (FWHM) sec/i-shaped pulse 
incidents into a uniform grating with 15.0 cm -1 wavenumber detuning from the center 
of the bandgap. The coupling strength of the fiber grating is 10 cm -1 , the nonlinear 
coefficient T is 0.018 cm/GW, and the group velocity v g is chosen to be c/n with n = 1.5 
and c being the speed of light in free space. When the input peak intensity is below the 
required value for forming a solitary pulse in the FBG (about 4.5 GW/cm 2 in this case), 
the peak intensity of the pulse will decrease along the propagation. On the other hand, 
when the input peak intensity is above 4.5 GW/cm 2 , the peak intensity of the pulse 
oscillates during the propagation within the grating. Only when the nonlinearity can 
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Figure 1. Dispersion relation of linear fiber Bragg gratings. 

exactly compensate the dispersion induced by the FBG, one can have a stable solitary 
pulse inside the grating. 

3. Quantum Fiber Bragg Grating Solitons 

After finding these classical solutions, we now turn to the calculation of the quantum 
properties. The Hamiltonian of the system can be expressed by the forward (U a ) and 
backward (U b ) field operators as follows [37] : 

H = - v g {6 J dz {UlU a + UlU b ) +ijdz (ul—u a - Ut-^Ub) (4) 

+ kJ dz (UlU b + UlU a ) 

+ \ J dz {UlUlUaUa + UlulUbUb) 

+ r Jdz {ululu b u a + ululuA)} 

Here the field operators U a and U b must obey the usual equal time Bosonic commutation 
relations. 

[U a (z 1 ,t),Ul(z 2 ,t)} = 5(z 1 -z 2 ) 
[U b ( Zl ,t),Ul(z 2 ,t)] =8{z x -z 2 ) 
[U a (z h t),U a (z2,t)} = [Ui(z u t),Ui(z 2 ,t)] = 
[U b ( Zl ,t),U b (z 2 ,t)} = [Ul( Zl ,t),Ul(z 2 ,t)] = 
[U a (zi,t),U b (z 2 ,t)} = [U a ( Zl ,t),Ul(z 2 ,t)} = 
In the Heisenberg picture we can derive the quantum nonlinear coupled mode 
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Figure 2. Measurement scheme of direct detection for observing FBG soliton 
squeezing. Here G is a gating device which will block out all the transmittive pulses 
but the first one; D is an optical detector. 



equations (QNCMEs) by the Hamiltonian in Eq. (J2J: 

--U a (z,t) + —U a = i5U a + mU b + iTU\U a U a + 2iTUllJ b U a (5) 

-Q- t U b {z, t) - —U b = i5U b + iKU a + iTUlU b Ub + 2iTU\U a U b (6) 

This derivation automatically proves that the QNCMEs preserve the commutation 
brackets. The form of QNCMEs in Eqs. (J5KJ) is the same as the classical NCMEs, 
except that it is now a set of quantum operator equations. For optical solitons with 
a large average photon number, we can safely use the linearization approximation to 
study their quantum effects. By setting 

U a (z, t) = U a0 (z, t) + u a (z, t) 
U b (z, t) = U b0 (z, t) + u b (z, t) 

where U a0 and U b0 are classical fiber Bragg grating soliton solutions of NCMEs, one 
can obtain a set of linear quantum operator equations that describe the evolution of 
the quantum fluctuations associated with the fiber Bragg grating solitons after the 
linearization approximation is performed. 

v a dt\u b [ iK + 2iVU*U w £ + 15 + 2iT\U a0 \ 2 + 2iT\U b0 \ 2 j 



+ 



iTU 2 aQ 2iTU aQ U w 
2iTU a0 U b0 iTU 2 




This set of equations is still a set of quantum operator equations for the quantum 
perturbation field operators u a (z,t) and u b (z,t), which have to satisfy the same equal 
time commutation relations as the unperturbed field operators U a (z,t) and U b (z,t). 

To solve the linear quantum operator equations in Eq. ((7|) by the back-propagation 
method, we define the inner product operation according to 



(/ID = \jdz irja + f a gl + fbdb + f b gl) (8) 

for a set of adjoint equations which satisfies the following property: 
<7 . . , >. 
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In other words, the inner product between the solutions of the two equation sets 
is preserved along the time axis. With this desired property, we can express 
the inner product of the output quantum perturbation operator with a projection 
function in terms of the input quantum field operators which have the know quantum 
characteristics. This will allow us to calculate the quantum fluctuations of any inner 
product between the output quantum operator with a given projection function. It is 
not difficult to show that the set of adjoint equations for the perturbed QNCMES in 
Eq. is given by: 

\ ( -§- z +i5 + 2iT\U a ,\ 2 + 2iT\U m \ 2 iK + 2iTU a0 U^ \, 

v g dt\u$ ) " \ in + 2iTU* a0 U b0 £ + iS + 2iT\U a0 \ 2 + 2iT\U b0 \ 2 ) \ 

( -iTU 2 -2iTU a0 U b0 
{ -2iTU a0 U b0 -iTU 2 

Actually under the linearization approximation the measurement of a physical 
quantity can always be expressed as an inner product between a measurement 
characteristic function and the perturbed quantum field operator [211! • By back- 
propagating the classical adjoint equation set with the measurement characteristic 
function as the initial condition, one can express the measured inner product as an inner 
product between the input quantum operator and the back-propagated measurement 
function. It is then quite easy to calculate the quantum uncertainties of the measured 
quantity if the quantum characteristics of the input quantum field operators are known. 
Typically we will assume the input quantum state to be a coherent state. The squeezing 
ratio of the measured quantity is defined according to: 

R(T) = var[(f(z)\u(z,t = T)] 
' ' var[(f(z)\$(z,t = Q)] 
var[(F T (z)\ii{z,t = 0)] 
~ var[(f{z)\u(z,t = Q)] 

Here var[-] means the variance, f(z) is the original projection function and Ft(z) is 
the back-propagated projection function. The choice of the characteristic function 
f(z) = (f a (z),f b (z)) T will depend on the measurement to be performed. When the 
measured quantity is the photon number of the pulse, the characteristic function f a (z) 
is simply the normalized classical output pulse and f b (z) = 0. If the homodyne 
detection is used, then f a {z) is the normalized local oscillator pulse. 

In Fig. |21 we illustrate the possible direct detection scheme for measuring the photon 
number squeezing of the FBG solitons. To avoid the complication due to the multiple 
transmitted pulses of the FBG, it may be necessary to use a time-gating device to make 
sure that only the first transmitted pulse is detected. 

Based on the formulation given above, we now apply the developed quantum theory 
to the case of a pulse incident into a uniform nonlinear FBG and calculate the photon 
number fluctuations of the first transmitted pulse. First of all, the transmittance of a 
50 cm long FBG versus different intensities of the input pulse is plotted in the top curve 
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Figure 3. Transmittancc (top) and photon number squeezing ratio (bottom) for fiber 
Bragg grating solitons with different input intensities. 

of Fig. El It can be seen that the transmittance of the FBG is a nonlinear oscillating 
function of the input pulse intensity. The calculated photon number squeezing ratio 
is shown in the bottom of Fig. El with the same parameters. One can see that the 
squeezing ratio decreases monotonically when the input intensity is below the intensity 
of the fundamental soliton. The optimum squeezing ratio occurs when the pulse energy 
of soliton is slightly large than that of the fundamental soliton. The curve begins to 
oscillate strongly when the intensity is much larger than that of the fundamental soliton. 
One important observation is that the oscillation behavior of the FBG transmittance and 
the squeezing ratio matches very well. That is, the squeezing ratio has a local minimum 
when the transmission has a local maximum. This agrees with the intuitive expectation 
that larger amplitude squeezing should occur when the transmittance curve is saturated. 
Another way to intuitively understand the results is as follows. The periodic grating 
structure acts like a spectral filter which can filter out the noisier high frequency 
components in the soliton spectrum and produce a net photon number (amplitude) 
squeezing effect just as in the previous soliton amplitude squeezing experiments where 
a spectral filter is cascaded after a nonlinear fiber I38j. 

To make sure that the FBG solitons are actually amplitude squeezed, we perform 
another calculation to simulate the squeezing ratio when the homodyne detection scheme 
is used and assume that the local oscillator pulse is exactly the classical output pulses. 
With the homodyne detection scheme, one has the additional degree of freedom to 
adjust the relative phase between the local oscillator and the signal. Fig. 0] plots the 
squeezing ratio for different input intensities and for different local oscillator phases 
when propagating through a constant FBG length (50cm 2 ). One can see that for small 
input intensities (below hGW/cm 2 ) the squeezing direction is close to but not exactly in 
the photon number (or amplitude) quadrature, 9 = 0. However, when the intensity of 
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Figure 4. Squeezing ratio of the quadrature field versus the input intensity and the 
local oscillator phase. 

the FBG soliton is large enough, the squeezing direction will be in the photon number 
quadrature. 

Fig. El shows the dependence of the photon number squeezing ratios for different 
FBG lengths with a constant input intensity (J = A.5GW/cm 2 ). We find that the 
squeezing ratio monotonically decreases with the increasing of the FBG length and 
saturates at the length around 60 cm. Intuitively the filtering effects of the FBG will 
unavoidably introduce additional noises on the light fields and will eventually cause the 
squeezing ratio to become saturated. In Fig. El we also plot the squeezing ratio for 
different FBG lengths and for different local oscillator phases with a constant input 
intensity (/ = 4.5GW/cm 2 ). Again when the FBG length is long enough, the squeezing 
direction will be in the in-phase (photon number) quadrature. 

4. Quantum FBG solitons in nonuniform gratings 

It is well known that one can engineer the dispersion along the FBG by using a 
nonuniform FBG. Such apodized nonlinear FBGs have been used for adiabatic soliton 
pulse compression within a very short length of several centimeters |34[ 135] . Intuitively 
the solitons with higher peak intensities will exhibit large squeezing due to higher 
nonlinear effects. It is thus expected that one should be possible to compress the 
pulsewidth of the FBG soliton and enhance its squeezing ratio simultaneously with 
the use of a suitably apodized FBG. To verify this idea, here we consider a nonuniform 
FBG which has a position dependent coupling coefficient described by 

k(z) = n + az (12) 

where kq = ujQe/2nc is the initial coupling coefficient and a is the slope of the coupling 
coefficient. 
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Figure 5. Photon number squeezing ratio for Bragg solitons after propagating through 
different length of FBGs. 




Figure 6. Squeezing ratio of the quadrature field versus the FBG length and the local 
oscillator phase. 

In the following calculation we consider the same 60 ps FWHM sec/i-shaped input 
pulse with the peak intensity of / = 4.5GW/cm 2 for the nonuniform grating without 
changing any parameter. In Fig. 0we plot the squeezing ratio versus the FBG length 
with a constant input intensity (J = A.5GW/cm 2 ) and different apodization slopes. 
The three curves in Fig. [7| correspond to different slopes of k(z) (a = +0.04,0, —0.04) 
respectively. When a = 0.04(l/cm~ 2 ), the FWHM of the original FBG soliton (60 
ps) will be adiabatically compressed down to 30 ps after propagating through a 70 cm 
FBG. Because of this, the achievable optimal squeezing ratio thus increases from 8 dB 
to 9.7 dB for the propagation distance of 70 cm. On the other hand, if the slope of 
the coupling coefficient is positive, we have a broaden FBG soliton and the optimal 
squeezing ratio will be degraded. If the slope of the coupling coefficient is too large, 
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Figure 7. Photon number squeezing ratio versus the FBG length with a constant 
input intensity (/ = A.hGW / cm 2 ). The coupling coefficient of the FBG is a function 
of propagating length: k(z) = kq + az. The three curves correspond to different slopes 
of k(z); Solid-line: a — 0; Dashed- line: a = —0.04; Dash-dotted line: a = +0.04. 

then the soliton cannot be compressed adiabatically and thus the optimal squeezing ratio 
will be degraded eventually. This can be seen in Fig. |HJ where we plot the squeezing 
ratio versus the grating length with larger slopes (a = —0.08, +0.08). The result for 
the uniform FBG case (a = 0) is also plotted for comparison. These results show that 
one can actually tailor the squeezing ratio of the FBG solitons as long as the soliton 
pulses evolve adiabatically. Recently some new theoretical results for the dynamics of 
classical gap solitons in an apodized grating have been carried out [36J. It was shown 
that the velocity of the gap soliton can be essentially slowed down in the apodized 
grating and interesting soliton collision behaviors including the merge of the solitons 
were predicted. We shall investigate more about the other quantum dynamics of FBG 
solitons in apodized fiber gratings in the future. 

5. Conclusion 

In conclusion, by following the line of theoretical development on the quantum theory 
of optical solitons initiated by Prof. Hermann A. Haus, we have developed a general 
quantum theory for bi-directional nonlinear optical pulse propagation problems and have 
especially used it to study the squeezing phenomena of fiber Bragg grating solitons. It 
has been shown for the first time that the output FBG soliton pulses will get amplitude 
squeezed automatically. The squeezing ratio of the FBG solitons exhibits interesting 
relation with the fiber grating length as well as with the intensity of the input pulse. 
The squeezing ratio saturates after a certain grating length and the optimal squeezing 
ratio occurs when the intensity of the FBG soliton is slightly large than that of the 
fundamental soliton. With the use of nonuniform FBGs, we also find that one can 
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Figure 8. Photon number squeezing ratio versus the FBG length with different slopes 
of k(z); Solid- line: a = 0; Dashcd-linc: a = —0.08; Dash-dotted line: a = +0.08. 
Other parameters are the same with those in Fig. 

compress the FBG solitons and enhance its squeezing ratio simultaneously, as long as 
the soliton pulses evolve adiabatically. To actually measure the quantum fluctuations of 
the fiber Bragg grating solitons experimentally, we propose to use a time-gating device 
to block out other smaller multiple transmitted pulses and only directly detect the first 
transmitted pulse from the grating. Since our calculation is already based on the existing 
experimental parameters, it shall be very interesting to see if one can actually verify the 
theoretical predictions experimentally. 
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